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Abstract: We construct a holographic model for the superconductor- insulator- 
superconductor (SIS) Josephson junction at zero temperature by considering a com- 
plex scalar field coupled with a Maxwell field in the four-dimensional anti-de Sitter 
soliton background. From the gravity side we reproduce the sine relation between 
the Josephson current and the phase difference across the junction. We also study 
the dependence of the maximal current on the dimension of the condensate operator 
and on the width of the junction, and obtain expected results. 
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1 Introduction 

With the help of the AdS/CFT correspondence [1], the superconductivity and super- 
fluid phenomena have been studied intensively on the gravity side in recent years [2- 
4]. For some reviews, see, for example, [5-7]. One of the interesting phenomena 
associated with superconductivity is the Josephson effect [8] (a brief description of 
the Josephson effect is given in appendix A), namely, the effect of electrons tunneling 
between two superconductors separated by a weak link. The weak link can be a nor- 
mal conductor or an insulating barrier. Correspondingly the Josephson junction is 
referred to as the superconductor-normal-superconductor (SNS) or superconductor- 
insulator-superconductor (SIS) junction, respectively. 

Recently, a holographic model for a three-dimensional SNS Josephson junction 
has been constructed and studied by Horowitz et al [9] in an AdS-Schwarzschild 
black hole background with a Maxwell field and a complex scalar field. The exten- 
sion to a four- dimensional Josephson junction has been discussed in [10, 11]. The 
Josephson junction array based on a designer multigravity has been constructed in 
[12], and the holographic p-wave Josephson junction has also been discussed in [13]. 
In those works, some familiar features of Josephson junction have been reproduced 
on the gravity side. Note that in those studies the dual weak link is a normal con- 
ductor. Therefore it would be of great interest to construct a holographic model for 
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a Josephson junction with an insulator link and study its features. This is the aim 
of the present paper. 

In order to build the holographic model for an SIS Josephson junction, we first 
need a model to realize the superconductor /insulator phase transition. Fortunately 
such a holographic model was constructed in [14]. There it was shown that in an 
Einstein-Maxwell-complex scalar field theory with a negative cosmological constant, 
one can have two phases. One is the AdS soliton solution with a vanishing scalar 
field. This solution is dual to a confined gauge theory with a mass gap in the AdS 
boundary, which resembles the insulator phase. The other is the AdS soliton solution 
with non-vanishing scalar field. When chemical potential increases beyond a certain 
value p c , the above AdS soliton solution turns out to be unstable, a new and stable 
AdS soliton solution with nontrivial scalar field appears. The new solution is dual to a 
superconducting phase. In this way one realizes the superconductor/insulator phase 
transition. The analytical study of the holographic insulator /superconductor phase 
transition in a five- dimensional AdS soliton spacetime was presented [15]. With the 
holographic model of superconductor/insulator phase transition, we here construct 
a holographic model dual to a (l+l)-dimensional SIS Josephson junction in a four- 
dimensional Einstein-Maxwell-complex scalar theory with a negative cosmological 
constant. 

In the probe limit, by numerically solving the coupled non-linear equations of 
motion for the Maxwell field and scalar field in an AdS soliton background, we show 
that the Josephson current is proportional to the sine of the phase difference across 
the SIS junction. In particular, it is found that the maximum current J max decreases 
when the mass square m 2 of the scalar field becomes large or the width L of the 
junction increases. In addition, the condensation (O) dual to the scalar field will 
also decrease exponentially with respect to L. The coherence length of the junction 
£ is found very close to each other from the fittings of the two fig UreS J max r ^ L 
and (O) ~ L. We also analytically study this holographic model by virtue of the 
Sturm-Liouville (S-L) eigenvalue problems [16]. We find that the critical chemical 
potentials p c obtained analytically are in good agreement with the ones from the 
numerical calculation. Besides, we analytically obtain the critical exponent 1/2 from 
the relation between the operator condensation and the chemical potential near the 
critical point, i.e., (O) oc yj ' p — p c ; The charge density p is found to be linearly 
proportional to the chemical potential, i.e., p oc (p — p c ), which is qualitatively 
consistent with the numerical calculation. 

This paper is organized as follows. In section 2, we set up the gravity dual 
model of an inhomogeneous superconductor in a four-dimensional Einstein-Maxwell- 
complex scalar theory with a negative cosmological constant. The holographic in- 
sulator/superconductor phase transition is studied in section 3. In section 4, we 
numerically solve a set of coupled non-linear equations of motion, and study the 
characteristics of the holographic SIS Josephson junction. Section 5 is devoted to 
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the conclusions. In appendix A, we briefly describe the Josephson junction, while in 
appendix B, we give an analytical study on the holographic superconductor /insulator 
phase transition. 



2 Setup 

Let us begin with the action of Einstein gravity with a cosmological constant in four 
dimensions: 

S= f d A x^(R -2A) , (2.1) 



where the negative cosmological constant A is related to I by A = —3/£ 2 , where £ 
is the radius of AdS space. In this system we have a Ricci flat AdS-Schwarzschild 
black hole solution as 

ds 2 = -^-dt 2 + £ 2 ^- + r\dx 2 + dy 2 ), (2.2) 
£ 2 f(r) 

where f(r)=r 2 — rQ/r and r = r$ is the black hole horizon. With the double Wick 
rotation, we can obtain the so-called AdS soliton solution of this system [17] 

ds i = _ r 2 dt 2 + f - T -pj 1 dr 2 + r 2 dx 2 + (V 2 - r -pj d X 2 ■ (2.3) 

Here the coordinate \ should have a period with (3 = Air£/3ro, otherwise there will 
be a conical singularity at r = tq. In that case r = tq in the metric (2.3) becomes 
a tip of a cigar-like geometry. The temperature associated with this AdS soliton 
vanishes. In addition, let us stress here that the coordinate x is infinitely extended, 
namely x G (— oo, oo). 

Next we consider the matter sector of the model: a Maxwell field coupled with 
a charged complex scalar field. Its action can be written as 



( --F^F^ - | - iA^\ 2 - m 2 \^\ 2 



where F^ u = d^A u — d u A^ is the Maxwell field strength and m is the mass of the 
scalar field. The equations of motion (EoMs) of the scalar and Maxwell fields read 

(V„ - iA a ) (V a - iA a ) i> - = , (2.5) 
V a F ah - i [^*(V b - iA h )i) - ^(V 6 + iA b )^*] = . (2.6) 

In order to solve the above equations, we choose an ansatz as 

^ = Me*\ A = (A t ,A r ,A x ,0) , (2.7) 

where ip, A t , A r , and A x are all real functions of coordinates r and x. Instead 
of A, we work with the gauge-invariant fields: M = A — d(p. 



- 3 - 



We are going to study the holographic model in the probe limit. That is, we can 
ignore the back reaction of the matter fields on the AdS soliton geometry (2.3). In 
that case, the EoMs for the matter fields in the AdS soliton background become 

1 / f 2\ 1 / M 2 M 2 \ 

#M + -r f dM + + m + j (-r - f M " ~ ~t - ™ 2 ) M = °> (2-8) 

1 2/ M \ ( f 2\ 

d r M r + —d x M x + — [M r dM\ + ^jdM) + {j + - r ) M <- = °> ( 2 - 9 ) 

1 f 2\ib\ 2 

d 2 r M t + -^dlMt + J jd r M t - -^-M t = 0, (2.10) 

d 2 x M r - d r d x M x - 2r 2 \ij\ 2 M r = 0, (2.11) 

d 2 r M x - d r d x M r + f j (d r M x - d x M r ) - 2 ^M X = 0, (2.12) 

where the superscript a prime denotes the derivative with respect to r. In the 
following, we will work with the case m 2 > —9/4, which is above the Breitenlohner- 
Freedman bound [18]. Here we have set £ = 1. 

Because Eqs. (2.8)-(2.12) are a set of non-linear coupled equations, one can- 
not solve these equations analytically. However, it is straightforward to solve them 
numerically. In order to numerically solve Eqs. (2.8)-(2.12), we need to specify the 
associated boundary conditions in the x and r directions. By imposing the Neumann- 
like boundary condition suggested in [14], we obtain the asymptotic behaviors of ip, 
M t , M r and M x near the tip r = ro as 

M -»• a (x) + ai (x)(r -r ) + O ((r - r ) 2 ) , (2.13) 

M t -)• 6 (x) + h(x)(r -r ) + O ((r - r ) 2 ) , (2.14) 

M r -> co(x) + ci(x)(r -r ) + O ((r - r ) 2 ) , (2.15) 

M x d (x) + rf x (x)(r - r ) + O ((r - r ) 2 ) . (2.16) 

Analyzing the EoMs of the scalar field and Maxwell field on the AdS boundary at 
r = oo, we can take the asymptotic forms as 

a ( ! \ (217) 

lrl r (3-v / 9+4m 2 )/2 r (3+V9+4m 2 )/2 \ r (3+v / 9+4m 2 )/2+l ' 



M r -> O ( 1 ) , (2.19) 



M B ->z/(x) + - + 0(^l , (2.20) 



where ip^ (i = 1,2) are the condensation values of the dual scalar operators (Oi), 
according to the AdS/CFT dictionary. We can set one of ip^ to be vanished. Here 
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we choose ip = 0. In addition, here p, p, v and J correspond to the chemical 
potential, charge density, superfluid velocity and current in the boundary field theory, 
respectively [19-24]. 

According to [14], in the above setup we can construct a holographic SIS Joseph- 
son junction by adjusting the boundary chemical potential p(x): in a region around 
x = along the direction x, p(x) is below the critical chemical potential p c so that 
this region behaves like an insulator; while beyond this region p(x) is above the crit- 
ical chemical potential so that those two regions outside the insulator region behave 
like two superconductors. In this way we construct a holographic SIS Josephson 
junction. Thus, the gauge invariant phase difference 7 across the insulator region 
can be defined as [9] 

/oo 
dx [v{x) - u{±oo)\ . (2.21) 
-00 

The Josephson current J across the Josephson junction has a sine relation to the 
phase difference 7. We will numerically confirm this in Sec. 4. 



3 Holographic superconductor/insulator phase transition 

In this section, we will first analyze the critical chemical potential p c for the su- 
perconductor/insulator phase transition near x — > 00, in order to give a proper 
configuration of the chemical potential for the SIS Josephson junction in the next 
section. In the limit x — > 00, we assume that all the fields are independent of x, i.e., 
all the fields are homogeneous in the x-direction. Therefore, Eqs. (2.8)-(2.12) can 
be reduced to two coupled ordinary differential equations for and as 



/' 2\ . ,. 1 (M x 



d 'M + + + j [it - ™ 2 ) w = , (3.1) 

f 2\ib\ 2 
d 2 r M t + J jd r M t - -^-M t = . (3.2) 

Here we will numerically solve the EoMs (3.1) and (3.2), while some analytical results 
will be presented in appendix B, for comparison. We will see that the critical chemical 
potentials obtained from the two approaches are very close to each other. 
The asymptotic forms of and M t near the boundary r = 00 are 

_ r (3-v / 9+4^)/i' + r (3+v / 9T4^)/2' ' ' ' 

M t = p- 9 -, (3.4) 
r 

in which the four parameters (V^ 1 ^, 4>^ 2 \ A*, p) are constants independent of x, p and 
p are the chemical potential and charge density, respectively, while ip^ and tp^ are 
the condensations of dual operators in the superconducting phase. 
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(a) (b) 

Figure 1. (a) The condensation of the operator (O) versus the chemical potential p; (b) 
The charge density p versus the chemical potential p. 

In the numerical calculations, it is convenient to use the shooting method to solve 
the boundary value problem for the ordinary differential equations (3.1) and (3.2). 
In the following we will set tq = 1 and £ = 1. In Fig. 1 we plot the condensation of 
the operator (O) and charge density p as functions of the chemical potential p. We 
can see that the charged scalar operator condensates when the chemical potential is 
above critical value p c : 

p c « 1.7182 for m 2 = -2, (3.5) 
p c « 2.2205 for m 2 = -5/4. (3.6) 

We can see that the critical chemical potential p c grows when m 2 increases, which 
means that the condensation becomes hard when m 2 increases. In the case of p < p c , 
the condensation does not happen and the charge density p is always vanishing, one 
realizes a holographic insulator phase; while p > p c , the charged operator conden- 
sates, which breaks the U(l) gauge symmetry spontaneously. Therefore, when the 
chemical potential crosses the critical value p c , one can realize a holographic insula- 
tor/superconductor phase transition. In appendix B, we get analytically the values 
of the critical chemical potential by using the S-L eigenvalue procedure, which are 
very close to the ones from the numerical calculations. 

In addition, near the critical point, we can read off the critical behavior of the 
condensation (O) and the charge density p from Fig. 1 as 

(O) « 1.3134 - 0c P ~ 1.4229(0 - p c ) for m 2 = -2, (3.7) 
(O) « 1.4231 y/p - p c , p « 1.2758(/j - p c ) for m 2 = -5/4. (3.8) 

We can see that (O) is proportional to the square root of p — p c , and p is linearly 
proportional to p — p c . We also obtain these square root and linear relations in the 
analytical study of the critical behavior in appendix B. 
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Figure 2. The sine relation between the Josephson current J and the phase difference 7 
for two different m 2 . Here the parameters are set to = 4, L = 3, e = 0.4, a = 0.5. 

4 Numerical solutions of the holographic SIS Josephson junc- 
tion 

In this section, we will numerically study the holographic SIS Josephson junction. 
For this aim, following [9], we choose the profile of the chemical potential fi(x) as 



/i(x) = jl c 



1 



2tanh(^ 



I j — tanh ( 



(4.1) 



where /ioo = yu(oo) = fi(-oo) is the value of the chemical potential \i at x — ±00, 
and the parameters L, a and e are the width, steepness and depth of the junction, 
respectively. Actually, we can suitably choose the parameters in fi(x) so that 

/jl(x) < yu c , for - L/2 < x < L/2, 

(4.2) 

fi(x) > /i c , for x < —L/2 and x > L/2. 

Thus the central part of the junction is a holographic insulator while the two sides 
are two holographic superconductors. Such a junction is an SIS Josephson junction 
as we expected. 

With the above chemical potential, we can numerically solve the set of coupled 
equations (2.8)-(2.12) by virtue of the spectral method. In Fig. 2, we plot the current 
J as a function of the phase difference 7 with /i^ = 4, L = 3, e = 0.4, a = 0.5. The 
circles and the squares are from the numerical calculations while the solid lines are 
the fittings of them. The relation between J and 7 can be fitted as 



J f=! 0.1336 SU17, for m 2 = -2, (4.3) 
J pa 0.0194 SU17, for m 2 = -5/4. (4.4) 
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(a) (b) 

Figure 3. (a) The maximal current J max versus the width L of the junction; (b) The 
condensation (0(0)) | j=o versus the width L of the junction, for two different m 2 . Here the 
parameters are chosen as = 4, e = 0.4, and a = 0.5. 



We see that the Josephson current J is indeed proportional to the sine of the phase 
difference 7, which realizes the SIS Josephson junction. Further we find that the am- 
plitude of the sine relation or the maximum current J max decrease when m 2 becomes 
large. The relation between the maximum current J max and the width of the junction 
L is plotted in Fig. 3(a); And the relation between the condensation (O(0)) = ip^ 
at zero current and L is given in Fig. 3(b). The circles and squares are obtained 
from the numerical calculations while the solid lines are the fittings of them. From 
appendix A, one learns that the dependencies of J max and (O(0))|j=o on L behave 
like 



L 



</ max = A e e , (4.5) 
(O(0))| J=0 = A ie -%, (4.6) 

where £ is the coherence length of the SIS Josephson junction, while A and A\ are 
two constants. In our holographic model, the fitted result in Fig. 3(a) gives 

J max w 1.4230e~ L/L2552 , for m 2 = -2, (4.7) 
J max w 0.7107e- L/a872 °, for m 2 = -5/4, (4.8) 

while in Fig. 3(b) the fitted result has the forms 

(O) w 2.1671e- L/(2xL4015) , for m 2 = -2, (4.9) 
(O) w 1.9783e- L/(2xa7969) , for m 2 = -5/4. (4.10) 

The coherence lengths £ fitted from the two plots are consistent to each other within 
10.5% error. Thus we see that in the holographic model, the familiar results (4.5) 
and (4.6) for the SIS Josephson junction can be indeed reproduced on the gravity 
side. 
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Figure 4. The maximum current J max versus /x c /Voo for two different m 2 . The parameters 
are set to L = 3, e = 0.4, a = 0.5. 

We also plot the relation between the maximum current J max and fi c / fi^ in 
Fig. 4, from which one can see that J max will decrease when \i c j //oo increases, and as 
He/ Poo goes close to 1, the current is suppressed very much. This result is expected, 
because when fi c /Hoc — > 1, the chemical potential in both sides of the insulator just 
reaches the critical value. In this case the superconducting current is very small and 
the "strength" of their superconductivity is rather weak. We expect that the current 
vanishes finally at fi c /Poo = 1- On the other hand, we can see the maximal current 
quickly increases when fi c / poo becomes small. This tendency is also understandable, 
since in this case the two sides of the junction are in the superconducting phases far 
from the critical point, and the "strength" of the superconductivity becomes very 
strong. 



5 Conclusions 

In this paper, we studied the holographic insulator/superconductor phase transition 
in a four- dimensional AdS soliton background. It was found that the critical chem- 
ical potential /i c increases when the mass square m? of the scalar field grows. We 
investigated the behaviors of the condensation (O) and the charge density p, and 
found that the numerical results near the critical point are in good agreement with 
those obtained from an analytical procedure. 

Furthermore, with the help of the holographic superconductor/insulator model, 
suitably setting the profile of the chemical potential fi(x) on the AdS boundary, and 
following [9], we constructed a holographic model for an SIS Josephson junction. We 
numerically solved a set of non-linear equations of motion and confirmed that the 
Josephson current J of the SIS junction is indeed proportional to the sine of the 
phase difference 7 across the junction. Besides, we also shown that the maximum 
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current J max decreases when the mass square of the scalar field becomes large. We 
further studied the dependence of J max and (C(0)) |j=o on the width L of the junction, 
and found that J max and (O(0))\j=q are suppressed exponentially when L increases. 
Thus we reproduced these main features of Josephson junction on the gravity side. 
In addition, J max was found to be a decreasing function of \i c j \i. 

The Josephson effect is an important phenomenon associated with superconduc- 
tor. Some features of Josephson junction are robust such as the sine relation be- 
tween the current and the phase difference across the junction. For the holographic 
models of the SNS Josephson junction constructed in [9] and the SIS Josephson 
junction discussed in this paper, one can indeed confirm this relation on the grav- 
ity side. Therefore, it is interesting to further study the holographic model for the 
superconductor-superconductor-superconductor junction, where the middle super- 
conductor is different from the other two superconductors. We expect that the sine 
relation will also hold in that case. Finally let us mention here that it is also of great 
interest to investigate the effect of magnetic field on the Josephson junction from the 
holographic point of view [33]. 
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A Brief description of Josephson junction 

To be complete, we will briefly review the Josephson junction in this appendix. Con- 
sidering two superconductors separated by an insulating barrier, in 1962, Joseph- 
son [8] made a remarkable prediction that there exists a current flowing across the 
middle insulator even when there is no voltage difference between these two super- 
conductors. This supercurrent is caused by quantum tunneling of Cooper pair. This 
is just the Josephson effect. The current is related to the difference of phases of 
macroscopic wave functions in the two superconductors: 

J = Jmax sin 7, (A.l) 

where J max is the maximal current (critical current) and 7 is the phase difference 
of the macroscopic wave function ip = \ijj\e ie in two superconductors, which is in- 
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troduced as an order parameter in the Ginzburg-Landau model. This remarkable 
prediction was experimentally confirmed soon [25]. 

The Josephson's prediction was based on a microscopic theoretical analysis of 
quantum mechanical tunneling of electrons through the insulating barrier layer. But 
in fact, this Josephson effect is much more general, it occurs when two superconduc- 
tors are connected by a weak link, the latter is called Josephson junction. The weak 
link can be an insulator as Josephson originally proposed, or a normal metal layer, 
or simply a short and narrow constriction. In those cases, the Josephson junction 
are called SIS, SNS and SCS junctions, respectively. 

Suppose that ipi = \tpi\e iei and ip2 — \ip2\e %d2 are two macroscopic wave functions 
describing these two superconductors. Due to the weak coupling between them, the 
two wave functions are not independent of each other. Instead they relate to each 
other as: d x ip\ = Kip2 and d x ip2 = —Kip\, where K is a parameter representing the 
strength of the weak coupling, and x is the direction perpendicular to the boundary 
surface between superconductor and insulator. In the case with two identical super- 
conductors, one has \ipi\ = j^l = = y/po, where po is the density of Cooper pair. 
Generally speaking, the supercurrent associated with a macroscopic wave function 
^(r) is given as J = — ij^r (\1/*V\[ / — ^V**), where e* and m* mean the charge 
and the mass of the Cooper pair. From this, one can obtain the Josephson current 
flowing over the weak link by substituting ipi = \ip\e iei and ^2 — | | e i6>2 into the 
supercurrent equation as 

Jx = J max sin(6 l 2 - 9i) = J max sin7, (A. 2) 

e*h 

along the direction x, where J max = K p is the maximal Josephson current. 

m* 

The maximal Josephson current J max is closely related to the width L of the 
weak link (see, for example, [26]) as 

Jmax ~ Po e" L/ «, (A.3) 

where the coherence length £, which is one of the characteristic scales of super- 
conductors, can vary with superconducting materials. The maximal current is also 
proportional to the density of Cooper pair in the weak link as J max ~ |^| 2 - Consider- 
ing the relation between the number density p and the condensation (O): p ~ (O) 2 , 
see, for example, (B.19), one has 

(O) - JUL - e' Lm , (A.4) 

in the weak link of the Josephson junction. Indeed, in our holographic model of SIS 
Josephson junction, we have reproduced the three relations (A.l), (A.3) and (A.4). 
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B Analytical study on the superconductor /insulator phase 
transition 



In this appendix we will take advantage of the S-L method [16] to analytically study 
the homogenous solutions of the equations of motion near the critical point of the 
phase transition. For further applications of the S-L method to other holographic 
superconductor models, see Refs. [15, 27-32]. 

B.l Analytical results of the critical chemical potential for general mass 

In this case all matter fields are independent of the coordinate x. By setting z = 1/r, 
the EoMs (2.8) and (2.10) are reduced to (we have set M t = <p here) 

*«*) - + = o. (B.i) 

- i^M') - ^rS)*( J ) = °- < B - 2 > 

At the critical point, the scalar field ip vanishes. Near the critical point, we can 
introduce a trial function F(z) into the asymptotical form of ip near the boundary 
like 

1>\»0~(O)z A F(z), (B.3) 

where A = 3/2 + a/9/4 + m 2 is the conformal dimension of the dual operator. It is 
easy to see from the boundary conditions of ip that F(0) = 1. For simplicity, we can 
set F '(0) = 0. Therefore, the simplest form of F(z) is F(z) = 1 — az 2 , where a is 
a constant. Besides, near the critical point of the phase transition, <fr(z) ~ fi. Thus, 
substituting formula (B.3) into Eq. (B.l), and multiplying K(z) = z~ 2{ - l ~ A \— l + z 3 ) 
to both sides, we can reach 



d z \K{z)d z F{z)j + - P{z) + Q(z)/2 2 j F(z) = 0, (B.4) 

with 

P(z) = -z~ 2 ^ 2 -^ [m 2 + A {3 - A(l - z 3 )}] and Q(z) = -z^ 1 -^. (B.5) 

Through the S-L eigenvalue method, we can obtain the critical chemical potential 
by minimizing the following functional at some value of a, 

f dz{K{z){d z F{z)) 2 + P{z)F{z) 2 } 

^ = k ji • (B.6) 

/ dzQ{z)F 2 {z) 
Jo 
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Figure 5. The critical chemical potential [Jiq clS cl function of the mass square m . The 
plot is obtained from minimizing the formula (B.6) for certain a. 



Fig. 5 plots the relation between \x c and m 2 . In particular, the critical chemical 
potentials for the cases of m 2 = —2 and —5/4 are, respectively, 



[i c « 1.71884 for m 2 = -2, 
fj, c & 2.22116 for m 2 = -5/4. 



(B.7) 
(B.8) 



We see that these values (B.7) and (B.8) of the critical chemical potential are in 
good agreement with the numerical results (3.5) and (3.6), respectively. 

B.2 Analytical relations of (O) ~ (fi — p c ) and p ~ (fj, — fj, c ) 

Near the critical point p — > fi c , the condensation value of the operator (O) is very 
small, therefore, we can expand <fi(z) in the series of (O) like 



cf>~V c +(0)((z), 



(B.9) 



where, ((z) is a coefficient function independent of (O). We can easily find that the 
behavior of <j)(z) near the tip of the background geometry imposes 



C(i) = o. 



(B.10) 



The EoMs of £ can be obtained from Eqs. (B.2), (B.3) and (B.9) by taking the 
leading order of (O) as 



Multiplying T(z) = — 1 + z 3 to both sides of (B.ll), we arrive at 



j- z (r{z)d z t{z] 



-2p c (0)z^- 2 F(z) 



(B.ll) 



(B.12) 
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Integrating both sides of the above equation, we can get the Neumann boundary 
condition of ( near z — > as 



(-l + z 3 )d z ((z)\l = d z ((z)\ 

= -2n c (0) f dz z 2A - 2 F(z) 
Jo 



( 1 2a a 2 \ .„ 

= - 2 ^ 0> (2A3T-2ATT + 2AT3j' < R13) 

Therefore, from Eq. (B.ll) and the boundary conditions (B.10) and (B.13) of (, we 
can obtain 

, 1 , n J l8a 2 z 2A+1 9z 2A 



4A 2 + 8A + 3 A-2A 2 



4A 2 + 8A + 3 V 3 3 

2Fi A?(A + i)i (2A+5);23 



2A 2 + 3A + 1 

18^ 2A+3 / 2A 2A , 

+ — : 2-Pi U, + 1; + 2; z 3 

4A 2 + 4A-3 I 3 3 



1 



6a 2 A (4A 2 - 1) ^ {0) 



2A 1 
~3~ + 3 



A(2A- 1)(2A+ l)(2A + 3) 
+2a 2 A(2A - 1) (y/3ir(2A + 1) + 9) - 4 v / 37raA (4A 2 + 4A - 3) 

-12aA (4A 2 + 4A - 3) ^ (0) (^p^) + (2^ A - 9) (4A 2 + 8A + 3) 
+6A (4A 2 + 8A + 3) ^ (0) ^ + 1^ + 3A (a 2 (4A 2 - l) + a (-8A 2 - 8A + 6) 
+4A 2 + 8A + 3^ (- log (z 2 + z + 1) + 2 log(l - z) - 2^ tan^ 1 P^t^ 



+2 7 + log(27) 



(B.14) 



where, 2-^1(0, 6; c, d) is the hypergeometric function, ip^ is the digamma function 
and 7 is the Euler gamma function. 

Near the boundary z = 0, we can expand <fr(z) from (2.18) and (B.9) as 

<j>(z) -v-pz-^+iO) U{0) + C'(0)z + ^"(O)^ 2 + ■•■)■ (B.15) 

Therefore, comparing the coefficients of the z° term of (B.15) we have 

/!-/i c ~(O>C(0)cx(O> 2 . (B.16) 
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From the critical values of p c and the corresponding a above, we can reach 



(O) « 1.08790V)" - He, for m 2 = -2, (B.17) 
(O) « 1.09588V/i - /i c , for m 2 = -5/4. (B.18) 

which are qualitatively consistent with the numerical results in (3.7) and (3.8). The 
critical exponent 1/2 between the condensation value and the chemical potential is 
also consistent with the one from the mean field theory. Besides, we notice from the 
z 1 terms of both sides of (B.15) that 

p ~ -(O)C'(O) oc (O) 2 oc(fi- fi c ). (B.19) 

This linear relation between the charge density and the chemical potential is also 
qualitatively consistent with the numerical results in (3.7) and (3.8). From (B.19) 
we obtain 

p « 1.05323(/i-/i c ), for m 2 = -2, (B.20) 
p « 0.95214(/i - fi c ), for m 2 = -5/4. (B.21) 

Here the discrepancy in the numerical factors from both methods in the relations 
(O) ~ y/fi — JT C and p ~ (fi — /i c ) might be due to the fact that we have only adopted 
a simplest form of the trial function F(z) in (B.3), if we introduce higher order terms 
of z into the trial function, this may improve the discrepancy of the numerical factors. 
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